Will black hole-neutron star binary inspirals tell us about the neutron star equation of state? 
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The strong tidal forces that arise during the last stages of the life of a black hole-neutron star binary may 
severely distort, and possibly disrupt, the star. Both phenomena will imprint signatures about the stellar struc- 
ture in the emitted gravitational radiation. The information from the disruption, however, is confined to very 
high frequencies, where detectors are not very sensitive. We thus assess whether the lack of tidal distortion 
corrections in data-analysis pipelines will affect the detection of the inspiral part of the signal and whether these 
may yield information on the equation of state of matter at nuclear densities. Using recent post-Newtonian 
expressions and realistic equations of state to model these scenarios, we find that point-particle templates are 
sufficient for the detection of black hole-neutron star inspiralling binaries, with a loss of signals below 1% for 
both second and third-generation detectors. Such detections may be able to constrain particularly stiff equations 
of state, but will be unable to reveal the presence of a neutron star with a soft equation of state. 
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Introduction. Coalescing binaries of stellar-mass compact ob- 
jects, i.e. black holes (BHs) and neutron stars (NSs), are a pri- 
mary target for gravitational wave (GW) searches performed 
with kilometer-size laser-interferometric detectors, such as 
LIGO and Virgo, and one of the most promising for a first 
detection. Binaries containing at least one NS are particularly 
interesting because their GW signal contains signatures of the 
physical conditions of matter at nuclear densities (see, e.g. Q]- 
9 1) and may thus eventually reveal the equation of state (EOS) 
of NSs, which is currently highly uncertain ifTOl . 

In the case of BH-NS binaries, in particular, the most rele- 
vant of such signatures was traditionally thought to be the GW 
frequency above which no hydrostatic equilibrium is possible, 
/tide- Because this frequency was identified with the com- 
plete tidal disruption of the NS, the GW amplitude was con- 
sequently assumed to decay rapidly for larger frequencies ifTTI 
and information on the EOS was thought to be inferrable from 
/tide f° r binaries in which the NS were to be disrupted be- 
fore the system reaches the Innermost Stable Circular Orbit 
(ISCO). This picture was considerably modified by the recent 
numerical -relativity simulations presented in [12| which, in- 
stead, highlighted that /tide effectively marks only the onset 
of mass-shedding. The complete tidal disruption, on the other 
hand, is achieved only later and at higher frequencies, so that 
the GW spectrum decays exponentially at a cutoff frequency 
/cutoff — (1-2 — l-5)/tide- Stated differently, in contrast to 
the expectations of 1111 . the GW spectrum shows no distinc- 
tive feature at the mass-shedding frequency /tide; conversely, 
information on the NS EOS is in principle accessible through 
/cutoff - The two frequencies do not coincide because the tidal 
disruption is not instantaneous and the GW frequency changes 
rapidly in time, so that the transition from a chirping signal to 
an exponentially decaying one is pushed to higher frequen- 
cies. The difference between /tide and / cu toff increases as 
/tide approaches the GW frequency at the ISCO, /isco US- 

If recent numerical relativity simulations have pointed out 
the importance of the cutoff frequency as a very significant 



marker of the NS properties, they have also highlighted that 
the accurate determination of such frequency may only be 
made through numerical simulations, the computational costs 
of which are still prohibitive for a complete scan of the space 
of parameters. Moreover, it will be the inspiral portion of the 
GW signal that will be fully contained in the most sensitive 
part of the sensitivity curve of ground-based detectors. 

In light of all this, two questions follow naturally. (1) Will 
the use of point-particle template banks for inspiral searches 
cause a loss in the number of detected signals because of the 
small but secular tidal effects that develop? (2) Will such ef- 
fects be large enough to provide us with information about the 
NS EOS? We address here both questions by computing the 
GW phase evolution of BH-NS binaries when tidal effects are 
taken into account and by contrasting it with the one expected 
when both the BH and NS are treated as point-particles. We 
find that the tidal corrections to the phase evolution depend 
sensitively on the NS EOS, but also that they are generally 
small. Hence, present data-analysis pipelines, which do not 
include tidal corrections, are sufficient for a successful de- 
tection of BH-NS inspirals, but only a particularly stiff EOS 
would mark the GW signal sufficiently for us to determine it. 

Methodology. Early inspiral tidal effects may be effectively 
described in terms of a single EOS-dependent tidal deforma- 
bility parameter A l5H7l [T3ll (also indicated as p,2 GJ0X de- 
fined as the ratio of the quadrupole deformation induced on 
a star and the static external quadrupolar tidal field inducing 
the deformation, in this case the tidal field of its companion. 
The tidal deformability parameter is calculated via a linear 
axisymmetric 1 — 2 static perturbation around the centre- 
of-mass of the star, along the axis connecting the two com- 
panions; it depends on the EOS via the NS radius, R, and 
the 1 = 2 apsidal constant, or Love number, hi, through 
the relation A = 2k 2 R 5 /(3G) EHEI- The £ = 2 tidal 
deformation is expected to be the dominant source of EOS- 
dependent modifications to the inspiral phase evolution up to 
the last few orbits of the inspiral prior to the merger, when 
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FIG. 1 . Tracks of inspiralling BH-NS binaries at 100 Mpc, for a — and a PS EOS (red lines), and their position with respect to the sensitivities 
of AdLIGO (black line) and AdVirgo (light blue line). To avoid cluttering, we show only the strongest and the weakest signals, which refer to 
a 1.9 Mq and a 1.2 Mq NS, respectively, in binaries with mass ratio q = 0.1 and q = 1/3. Each track is terminated at (/end, h(f w i) V7eo&)> 
the red dashed continuations are shown only for reference. The shaded region is the one spanned by all the termination points obtained by 
varying Mns and q. A magnification reporting the lines at constant Mns and q that limit the shaded area is shown in the right panel. 



the mass-shedding takes over; higher-order in £ and A, non- 
linear, and viscous dissipation corrections are all considerably 
smaller [6]. The imprint of A on the inspiral of binary NSs 
was investigated in [6|, which concluded that, at a distance of 
100 Mpc and using the inspiral below 450 Hz, AdLIGO would 
be only be able to constrain A for an extremely stiff EOS. 

In order to estimate the GW phase accumulated during the 
inspiral of BH-NS binaries due to tidal effects, we follow |6 
[T31 . with the obvious difference that only one NS provides the 
tidal contribution. We also include the first post-Newtonian 
(1PN) order tidal corrections calculated in ifTTl . Finally, while 
in (6) the phase evolution was truncated at 450 Hz (roughly 
80% of /isco for two NSs with M NS = 1.4 M and R = 
15 km), we here set the value of the truncation frequency 
/end to be the smaller of /tide and /isco- The latter is the 
frequency yielded by the minimum of the PN energy func- 
tion (see below), while f t ^ e is determined using the relativis- 
tic toy model recently discussed in fl8l . which reproduces 
numerical -relativity results for fade within < 1%. In such 
a model the NS is described as an ellipsoid which deforms 
during the inspiral under the effect of the BH tidal field and 
its internal forces; this is consistent with the approach of (6), 
since only the £ = 2 tidal deformations are considered. 

We contrast the point-particle description of the quadrupole 
GW phase (<fipp) with the one obtained for a deformable, 
finite-size NS (<fi\) by integrating the tidal corrections to 
d<p/dx, where x is the PN expansion parameter, up to 1PN 
(relative) order IfTTl . i.e. we compute A<f) = <fipp — 4>\ (which 
is positive, due to the attractive nature of the tidal coupling). 
We moreover determine the quadrupole mode of the gravita- 
tional radiation in the frequency domain h(f). We use the 
inspiral part of spinning binary hybrid waveforms [Sect. Ill 
of |l9l l, augmenting the TaylorF2 phase ijj with A-dependent 



tidal corrections up to 1PN (relative) order [Eq. (3.9) in flD l. 
The PN expressions in |[T9l are derived from a 3PN expan- 
sion of the binding energy of the system and a 3.5PN (rela- 
tive order) expansion of the GW energy flux [Eqs. (3.3)-(3.4) 
in lfl9l l. within the commonly used stationary phase approx- 
imation [20]; spin terms are included up to 2.5PN order in 
phase and 2PN order in amplitude. We use the TaylorT4 ETI 
Description to replace the dx jdt in the expression of the am- 
plitude A. Even though leading-order tidal contributions are 
of 5PN order, as opposed to the point-particle terms which are 
kept up to 3PN or 3.5PN order, their coefficients "compen- 
sate" the smallness due to their PN order, making them com- 
parable to 3PN and 3.5PN terms (see lTl4l and Q for detailed 
discussions). We perform all integrations in the frequency do- 
main, from / s t alt = 10 Hz (roughly the low-frequency cutoff 
of AdLIGO/AdVirgo) to / end — niin(/ti(je, /isco)- 

The free parameters of our model are the NS (barotropic) 
EOS, the NS mass Mns, which ranges between 1.2M and 
the EOS-dependent maximum value A/ max , the binary mass 
ratio q = M N s/AI Bll , which we vary between 1/10 and 1/3, 
and the dimensionless BH spin parameter a, which is taken to 
be (anti-)aligned with the orbital angular momentum (equato- 
rial inspirals). To bracket the possible impact of the EOS on 
BH-NS inspirals and their GW radiation, we consider three 
EOSs dubbed APR, GNH3, and PS, respectively. Our choice 
is motivated by the fact that the APR EOS 1 22 1 is based on nu- 
clear many-body calculations, may be favoured by the obser- 
vations 121], and yields NSs with low A; the GNH3 EOS l24l . 
on the other hand, is based on mean-field theory and yields in- 
termediate A; finally, the "liquid" version of the PS EOS E5ll 
is somewhat dated but yields rather high A. The set chosen, 
therefore, covers the relevant range of stiffness (APR being 
the softest and PS the stiffest) and can be considered represen- 
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FIG. 2. Tidal distortion contribution to the quadrupole GW phase 
A(f> for the three representative EOSs. The tracks end at / en d and 
yield larger dephasings for stiffer EOSs, such as PS. 

tative of a much larger sample of EOSs. Note that all of these 
EOSs have a maximum mass 1.93 < M max / Mq < 2.66. 

Figure [T] illustrates the role played by / en( j. We consider 
the PS EOS and non-spinning BHs and show the tracks of in- 
spiralling BH-NS binaries at 100 Mpc (red solid lines), along 
with the sensitivity curves of AdLIGO (black line) and Ad- 
Virgo (light-blue line). The signal amplitudes are averaged 
over sky location and relative inclination of the binary. We 
show explicitly only the strongest and weakest signal, which 
refer to M NS = 1.9 M , q = 0.1 and to M NS = 1.2 M , q = 
1/3, respectively. The tracks terminate at (/end, M/end) V/end) 
and their continuations as red dashed lines serve only as a ref- 
erence. The shaded region, which is magnified in the right 
panel, is the one spanned by the termination point for all com- 
binations of M NS /M e [1.2, 1.9] and q e [0.1, 1/3]. 
Dephasing and Overlaps. Once a binary with parameters 
(q, a, -Mns, A) is selected, we compute A0(/ en d). Overall, we 
find that A<fi: (i) is greater for bigger A's, i.e. for more de- 
formable NSs (see Fig. [2); (ii) grows with q, i.e. for com- 
parable masses; (Hi) decreases as Mns [cf- Eq. (21) in (6)]; 
(iv) depends only weakly on the BH spin, since the only 
spin dependence may come through / en d, but binaries with 
/end = /tide < /isco are hardly affected, since / tid e is not very 
sensitive to a, while binaries with / en d — /isco < /tide are 
those with high Mns and low A, so that the gain or loss in 
/isco does not modify A</> significantly. 

To determine whether the dephasings found may affect the 
detection of BH-NS inspiral events, we compute, for each 
binary, the overlap between the point-particle model of the 
GW inspiral signal (/ipp) and the one which includes tidal 
deformability effects (h\); this is the normalized inner prod- 
uct of the two signals, maximised over time and phase shifts, 
i.e. 0[h pp , h x ] ee max {Wo} ^^h^) ■ where 

inner product is (h pp \h x ) = 4jft df ^s^tf} W . S h (f) 
being the noise power spectral density of a chosen detector. 
Note that we are implicitly assuming that the waveforms in- 
cluding tidal effects are the "real" signals and treating the 
point-particle waveforms as the templates used to detect them. 

Our results for the three EOSs considered and a BH with 
a = are shown in Fig.[3]for AdLIGO. Note that for any EOS 
choice and for any combination of the BH and NS masses, the 
overlap is always greater than 0.997, which corresponds to a 
1% loss of signals; this is true even for spins up to a = 1. 
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FIG. 3. Overlaps between PN waveforms for BH-NS binary systems 
modelled as point-particles ("PP") and with the inclusion of tidal dis- 
tortion effects ("A"). The overlap is calculated for the AdLIGO de- 
tector and for non-spinning BHs. 

The smallest overlap is given by the PS EOS combined with 
Mns = 1-2 Mq, q — 1/3, and a = 1 (inclusion of spin 
changes overlaps by < 10~ 3 ). Hence, even if all binaries were 
to have these extreme properties, the loss of signals would be 
less than 1%. All in all, BH-BH inspiral templates will allow 
second-generation interferometers to detect inspiralling BH- 
NS binaries with less than a 1% loss of signals. Similar re- 
sults hold for the third-generation detector Einstein Telescope 
(ET) |26|, with a minimum possible overlap of about 0.995. 
Measurability. Determining that the fraction of lost signals 
is below 1% does not address the question of whether the 
detected signals may be used to learn about the EOS. To 
address this point, we consider a nominal detector-binary 
distance of 100 Mpc and calculate the distinguishability as 
Shpp x = (hpp — h\\h ?P — h\) > 2(1 — 0[hp P , h\])p 2 , where 
p 2 = (hpp\hp P ) ~ (h\\h\) is the signal-to-noise-ratio (SNR) 
and we neglected the term (/ipp|/ipp) — (h\\h\) ^10~ 4 . Since 
we treat hpp as the "template" and h\ as the "signal", a neces- 
sary (but not sufficient) condition to distinguish between the 
two is that Shpp \ > 1 l2~71 . Clearly, the greater 5h PP \, the 
higher the chances of measuring the tidal effects. 

In Fig. [4] we consider AdLIGO and report 5hpp.\ for bi- 
naries with a = 0. In calculating the distinguishability, we 
average the signal amplitudes both over sky location and over 
relative inclination of the binary. Note the existence of a re- 
gion where Sh PP \ < 1 for all EOSs (white area) in the space 
of parameters (q, Mns); this indicates that the inspiral of BH- 
NS binaries falling in this region will not be distinguishable 
from a BH-BH inspiral. For larger mass ratios and smaller NS 
masses, Slipp y increases, becoming equal to 1 first for the PS 
EOS (red-shaded area) and then for the less stiff GNH3 EOS 
(blue-shaded area). The maximum value of 5hpp,\ is ~ 10 
(~ 5) for the PS (GNH3) EOS. Note that the black star pin- 



4 



1.9 




FIG. 4. Distinguishability <5/ipp,a by AdLIGO for binaries at 
100 Mpc and with a = 0. In the white area 5/ipp,A < 1 for all 
EOSs, in the red-shaded one <5/tpp,ps > 1, and in the blue-shaded one 
<5/f-pp,GNH3 > 1. Contours for /ipp,a = 1, 4 are shown for both EOSs. 
The black star marks a canonical 10 Mq-1.4 Mq BH-NS binary. 

points a canonical 10 Af Q -1.4 Mq BH-NS binary and that, for 
the popular APR EOS, (5/ipp,APR < 1 f° r an y (<Z, A/ns)- How- 
ever, an APR(-like) EOS may become measurable if the bi- 
nary is optimally oriented and/or is less than 100 Mpc away. 
The PS, GNH3, and APR EOS are (marginally) distinguish- 
able by AdLIGO within a range of 310, 225, and 75 Mpc, re- 
spectively, for a 3.6 M -1.2 M BH-NS binary. 

To test the robustness of our results, we repeated all calcu- 
lations after artificially decreasing / en( j by 20%. Even though 
the overlaps (distinguishabilities) increased (decreased) a lit- 
tle (e.g. the minimum overlap increased by less than 10~ 3 and 
the regions in which h PP ^ > 1 were practically unmodified), 
our conclusions remained unchanged, thus suggesting that our 
results are robust even for non-negligible changes of / en( j. 
Conclusions. The use of point-particle templates to detect BH- 
NS inspirals leads to a loss of detected signals which is below 
1% for both AdLIGO/Ad Virgo and ET. Moreover, for bina- 
ries at 100 Mpc, AdLIGO/ AdVirgo will essentially be "blind" 
to tidal effects if superdense matter follows an APR-like EOS 
and may be able to reveal them only for NSs with a partic- 
ularly stiff EOS and in large-mass-ratio binaries. This sce- 
nario improves for a more sensitive detector such as ET. In 
this case, the larger SNRs lead to a two-order of magnitude 
gain in 5hpp t \, so that even soft EOSs, like APR, are distin- 
guishable. 

A final remark should be made. The results presented here 
make use of the most accurate PN expressions available to 
date and the comparisons made with and without tidal correc- 
tions remove in part the problem of systematic biases in the 
PN formulation, which could be large for q < 1. Yet, they do 
not account for higher-order deformation corrections which 
could amplify tidal parameters [7 , 9|, and nonlinear responses 
to the tidal field, such as those produced by crust fracturing 



or resonant tidal excitation of stellar modes [28 1. Although 
these contributions should yield only fractional corrections to 
the already-small dephasings, either increasing or decreasing 
them, it is only their proper inclusion that will provide con- 
clusive limits on the tidal effects that affect the inspiral signal. 
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